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KUMARASWAMY LINDLEY-POISSON DISTRIBUTION: THEORY 
AND APPLICATIONS 


MAVIS PARARAI, BRODERICK O. OLUYEDE AND GAYAN WARAHENA-LIYANAGE 


ABSTRACT. The Kumaraswamy Lindley-Poisson (KLP) distribution which is an ex- 
tension of the Lindley-Poisson Distribution [21] is introduced and its properties are 
explored. This new distribution represents a more flexible model for the lifetime data. 
Some statistical properties of the proposed distribution including the shapes of the 
density and hazard rate functions are explored. Moments, entropy measures and the 
distribution of the order statistics are given. The maximum likelihood estimation 
technique is used to estimate the model parameters and a simulation study is con- 
ducted to investigate the performance of the maximum likelihood estimates. Finally 
some applications of the model with real data sets are presented to illustrate the 
usefulness of the proposed distribution. 


1. INTRODUCTION 


Lindley distribution [16], studied by Lindley in the context of fiducial and Bayesian 
statistics, is very useful for modeling failure time data. This distribution accommodates 
hazard rate functions that are increasing, decreasing or constant. However, models with 
complex hazard rate shapes such as unimodal, bathtub and other shapes are desirable 
in reliability analysis, human mortality studies and related areas. Ghitany et al. [11] 
proposed and studied the power Lindley distribution. Properties and applications of 
the Lindley distribution have been studied in the context of reliability analysis by 
several authors including Ghitany et al. [10], Sankaran [26] and Asgharzadeh et al. 
[1]. Nadarajah et al. [20] proposed and developed the mathematical properties of 
the generalized Lindley distribution. Properties of the exponentiated Power Lindley 
distribution were studied by Warahena-Liyanage and Pararai [35]. 

Several new families of distributions have been derived by compounding the Poisson 
distribution with many other continuous distributions to provide more flexible distribu- 
tions for modeling lifetime data. Kus [15] studied the exponential-Poisson distribution. 
Lu and Shi [17] derived and studied the Weibull-Poisson distribution. The exponen- 
tiated Weibull-Poisson distribution which generalizes the Weibull-Poisson was studied 
by Mahmoudi and Sepahdar [18], and the beta Weibull-Poisson was introduced and 
studied by Percontini et al.[23]. Barreto-Souza and Cribari-Neto [3] studied the ex- 
ponentiated exponential-Poisson distribution. The two parameter Poisson-exponential 
distribution with an increasing failure rate was studied by Cancho et al.[5]. Recently, 
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Pararai et al. [21] studied the properties of the exponentiated power Lindley-Poisson 
distribution, thereby generalizing the Lindley-Poisson distribution. 

The properties of Kumaraswamy [14] distribution were explored in detail by Jones 
[13]. The author contrasted the Kumaraswamy distribution with the beta distribution. 
Some of the good properties of the Kumaraswamy distribution include a simple nor- 
malizing constant, closed form solutions of the distribution and quantile functions as 
well as simple formulas for the moments. Cordeiro et al. [7] studied the Kumaraswamy- 
Weibull distribution and applied the model to some failure data. 

Motivated by the advantages of the generalized distribution with respect to having a 
hazard function that exhibits increasing, decreasing and bathtub shapes, as well as the 
versatility and flexibility of compounding Lindley and Poisson distributions in model- 
ing lifetime data, we propose and study a new distribution called the Kumaraswamy 
Lindley-Poisson (KLP) distribution, which inherits these desirable properties that also 
cover the shapes of quite a large number of models. 

We are also motivated to study the KLP distribution because of the wide and ex- 
tensive usage of Lindley distribution and the fact that the current generalization still 
provides a useful means for its continuous extension to more complex situations. An 
important and positive point of the current generalization is that the Lindley distribu- 
tion is a basic model or exemplar of the proposed KLP distribution. 

This paper is organized as follows. In section 2, the Kumaraswamy-G distribution 
(see Cordeiro and de Castro[8] for additional details), the model, its sub-models and 
some statistical properties including expansion of density function, quantile function, 
hazard function are presented. In section 3, we present the moments. Section 4 contains 
the distribution of the order statistics and Rényi entropy. Mean deviations, Bonfer- 
onni and Lorenz curves are presented in section 5. Maximum likelihood estimates of 
the model parameters and asymptotic confidence intervals are given in section 6. A 
simulation study is also presented in section 6. Section 7 contains applications of the 
proposed model to real data, followed by concluding remarks in section 8. 


2. THE MODEL, SUB-MODELS AND PROPERTIES 


The probability density function (pdf) and the corresponding cumulative distribution 
function (cdf) of the one-parameter Lindley distribution [16] are given by 








2 
(2.1) f(x; 6) = £ 7 +aje**, 2>0,6>0, 
and 
(2.2) Fis) =1— (: re ) eB 
Behl ‘ 


for x > 0,a, 6 > 0, respectively. 

Suppose that the random variable X has the Lindley distribution where its pdf and 
cdf are given in equations (2.1) and (2.2). Given N, let Xj,...,X~ be independent 
and identically distributed random variables from Lindley distribution. Let N be 
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distributed according to the zero truncated Poisson distribution [6] with pdf 
ore? 

TE 

Let X=max(Y},..., Yy), then the cdf of X|N = n is given by 


PIN =n) = [eal re’ emera oa 


GxiN=n(£) = h — (1 + F) aul , x>0,8>0,0>0, 


which is the exponentiated Lindley distribution. The Lindley-Poisson (LP) distribution 
denoted by LP(8,0) is defined by the marginal cdf of X, that is, 


Bax —pr 
t-ep foli- (1+ #5)e f |} 
(2.3) Grrl 6,8) a 

—e 
for x > 0, 8 > 0,0 > 0. The LP density function is given by 


(2.4) giP(x; 8,0) = sale a ua i fi) ai f 








for x > 0,8 >0,0 > 0; 


2.1. Kumaraswamy Lindley Poisson Distribution. In this sub-section, we present 
the Kumaraswamy Lindley-Poisson (KLP) distribution and derive some of its proper- 
ties including the cdf, pdf, expansion of the density, hazard function, quantile function 
and sub-models. 

Consider G(x) to be an arbitrary baseline cdf in the interval (0,1). The cdf G(x), 
referred to as Kumaraswamy-G distribution [8] has cdf 


F(z;a,b) = 1-(1-G(z)%?, 


where a and b are shape parameters. The pdf of the Kumaraswamy-G distribution is 
given by 


(2.5) f(aj;a,6) = abg(x)[G(z)]* [1 — Giz)? 1, a>0,b>0, 


where g(x) = a is the pdf corresponding to the baseline cdf. 


By taking G(x) as the cdf of the Lindley-Poisson (LP) distribution in equation (2.5), 
we obtain the Kumaraswamy Lindley-Poisson (KLP) distribution with a broad class of 
distributions that may be applicable in a wide range of day to day situations including 
applications in medicine, reliability and ecology. The cdf of the four-parameter KLP 
distribution is given by 


Os Harasani -efet (1+ ere] ii 


1— e? f 
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for x >0,0>0,8>0,a>0,b>0. The corresponding KLP pdf is given by 
2 —Ba = Bx \ e-r 
ab0 6? (1 + x)e™”? exp foj: (1+ a) |} 











fxrp(@) = cera 
1 -ep {61 _ (1+ aeh a—1 
x ea 
ao heete] 


for x > 0,8 > 0,0 >0,a >0,b >Q. 

Plots of the pdf of KLP distribution for several values of 6,0,a and b are given in 
Figures 2.1 and 2.2, respectively. The plots show that the KLP distribution is right 
skewed and can be decreasing (L shaped). 
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FIGURE 2.1. Plot of the PDF for different values of 6,0,a and b 
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FIGURE 2.2. Plot of the PDF for different values of 6,0,a and b 


Bx 
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Using the substitution ÀA = A(x) = 0 i — (1 + 
of the KLP distribution as 

f= ab B?2(1 + zee (1—e\*" i E a a lela 

KIPO (B+ 1)(e? -—1) \1-e? 1— e? 


2.2. Expansion of the Density Function. The expansion of the pdf of KLP distri- 
bution is presented in this sub-section. For b > 0 a real non-integer, we use the series 
representation 


Pi we can write the pdf 





(1 = [Grp(2)|*)?* = D (’ 7 ') (-1) (Grp), 


j=0 \ J 


Grela) = Grele: 8,8) = — ul 7 C y cea aa = 


e? —1 


where 
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We can rewrite the density of the KLP distribution as 


(oe) 


fxip(t) = (1p i p ') abgrp(x)|Gip(a)]7re* 


j=0 


ab08?(1 + x)e~** exp fo h — ( + fr) aa \ 
(8 + 1)(e® — 1) 




















7 aj+a—1 
= 1 — exp {6|1~ (1+ g) 
/b—1 B+1 
x Den(?5") — 
j=0 
oo fad eee ab) Ca) [e%(*+1) — 1] 
> -DEF 
B?(1 + x)e™®" exp {oth +1) h — (1 + A) \ 
á (Pea 
= p Wy ng (X; B,A(k +1)), 
j,k=0 
where 
be aende a 
. DR TEA a Ptah) a 


(e? —1)(k +1) 

and g(x; 6,0(k+1)) is the Lindley-Poisson pdf with parameters 6 > 0 and 6(k+1) > 0. 
This shows that the KLP distribution can be written as a linear combination of Lindley- 
Poisson density functions. Hence mathematical properties of the KLP distribution can 
be obtained from those of the LP properties. 


2.3. Survival and Hazard Rate Functions. The hazard function for the KLP dis- 
tribution will be presented in this sub-section. Using some selected values of 8,0,a, 
and b, some plots of the hazard function are presented. The hazard function of the 
KLP are given respectively by 


fxip(x; p, 0, a,b) 











hrel) 7 1 — Fgrp(z; 6,0, a,b) 
ab aere eee E 1- e\ 7! 
(6 + 1)(e? — 1) 1— e? 1— e? 
for x > 0,8 > 0,8 > 0,a > 0 and b > 0, where A = 6} 1 — (1+ Jem. The 


graph of the hazard function for various values of the parameters 8, 0,a and b are given 
in Figures 2.3 and 2.4, respectively. These graphs show that the KLP distribution is 
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suitable for monotonic and non-monotonic hazard behaviors which are more likely to 
be encountered in real life situations. 
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FIGURE 2.3. Plot of the Hazard Function for different values of 6, 6,a 
and b 
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FIGURE 2.4. Plot of the Hazard Function for different values of 8,0,a 
and b 


The graph of the hazard function for different values of the parameters exhibits 
various shapes such as monotonically increasing and bathtub shapes. 


2.4. Some sub-models of the KLP distribution. In this sub-section, some sub- 


models of the KLP distribution for selected values of the parameters (,6,a and b are 
presented. 


e When a =b = 1, we obtain the Lindley Poisson (LP) distribution whose cdf 
and pdf are given in equations (2.3) and (2.4). 


e When b = 1, we obtain the exponentiated Lindley Poisson (ELP) distribution 
which belongs to the resilience parameter family. 


e When a = b = 1 and 0 > 01, the Lindley, (L) distribution becomes a limiting 
form of the KLP distribution. 
e When b = 1 and 6 > 0°, we get the exponentiated Lindley (EL) distribution. 


e When a = 1, we get the new lifetime distribution belonging to the frailty 
parameter family 
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2.5. Quantile Function. The quantile function of the KLP distribution is obtained 
by solving the equation F'(Q(p)) = p, where 0 < p < 1. We therefore have 


aq b 
1 — exp {6 h — (1 + sa ) ezra] \ 
1— b= = 


1— e? 








alie exp fofi + a Í E 
so that 
zoezo = SF -znj a- Aa- a-p], 
thus, 
Z(p) =W (ah fı — sin h —~(1-—e*)(1-(1 = py] \) 


for 0 < p < 1, where W(.) is the Lambert W function [9]. The quantile function of the 
KLP distribution is obtained by solving for Q(p) in the above equation to obtain 


an = ai- pe eee) 
(2.9) 





3. MOMENTS 


In this section, we present the moments of the KLP distribution. Moments are 
necessary and important in any statistical analysis, especially in applications. They 
can be used to study the most important features and characteristics of a distribution 
(e.g., tendency, dispersion, skewness and kurtosis). 

The rt” moment of a random variable X following the KLP distribution, denoted by 


Lt, is 





wi kO B? m r —pr 
= a (1+ x)e 
XL PFDA 7 J, 


(3.1) x exp {otk +1) i = (: + a)l har. 


where wj% is defined in equation (2.8). In order to find the moments, consider the 
following lemma: 
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Lemma 1. Let 


Lı(8,0(k + 1), a,b, r) = [ a”(1+ 2)e** exp {otk +1) h = (1 + F) a hac 








ar 
then 
o ām p qtl 
m\ (p\ /q + 1\ (-1)?0™(k + 1)™81 
Li(8,0(k + 1),a,b,r) = ( a ) 
2 2 g0 s0 NPZ AN S m!(8 +1) 
. T(r+s+4+1) 
[B(p + 1)]rtstt 
Proof. Using the series expansion, e? = > wm we can rewrite the above integral as 
p=0 


I1(6,0(k4+1),a,b,r) = DA [vate fi (14 SE ee “| ae 


- EEZE OO) Sar 





q 











By letting u = (p + 1)x, we have x = Api and dr = Thus 
oe I(ir+s : 
ah tse PrP) gg = : 
/ [G(p + ])|erett 














By using Lemma 1, the rt” moment of the KLP distribution is 


0 
aD H T bn) 


j=0 k=0 





4. ORDER STATISTICS AND RÉNYI ENTROPY 


In this section, the distribution of order statistics and Rényi entropy for the KLP 
distribution are presented. The concept of entropy plays a vital role in information the- 
ory. The entropy of a random variable is defined in terms of its probability distribution 
and can be shown to be a good measure of randomness or uncertainty. 
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4.1. Order Statistics and Entropy. Suppose that X1,--- , Xn is a random sample 
of size n from a continuous pdf, f(x). Let Xin < Xan < +++ < Xn denote the cor- 
responding order statistics. If X1,--- , Xn is a random sample from KLP distribution, 
it follows from equations (2.6) and (2.7) that the pdf of the kt? order statistic, say 
Yk = Xk:n is given by 

n-k 


n! fietp(«) 
falus) = ED Hl 





PDFer 


Aa 
saa 2) > r 


j=0 m=0 p,q=0 


i a ') (— ])itm+pt+qtapta— lable? (q+1) _ 1] 
q (FP — 1)? (G41) 
O(q + 1)67(1 + we enr 

(B+ 1)[e%q+)) — 1] 


bp Orpen as 0 Gs b)g(x; 6, 0(q + 1)), 


= IÈ ie ogee | ar ea 


Me m,p=0 q=0 
ap +a-— 1 ae = 1] 
wt q 5 (e? — Leen (gl) i 











where 


Pim pal b, 0, a, b) 








and g(x; 3,0(q4+1)) is the LP pdf with parameters 8 > 0 and 6(q¢+1)) > 0. Thus, the 
distribution of the kt? order statistic is a linear combination of the LP distribution. 


4.2. Rényi Entropy. Rényi entropy [24][25] is an extension of Shannon entropy [29][30]. 


Rényi entropy is defined to be A, (fxrip(x; 8,0,a,b)) = log(Jo EN where 
v > 0, and v # 1. Rényi entropy tends to Shannon entropy as v z 1. We iera 
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have 


Ay(fxie) = — oe Gs firele) )]| 


s N Clay 


i, j,k=0 m=0 p=0 q=0 
(—1)itwtai-vtjtm o (j + v)]" Zf ate Ped 
0 

















i 2 FMR eo — 1rtai—v 

= a a 
OOC eee 

i (z sie 5) FE D] l 


for v > 0,v Æ 1. 


5. MEAN DEVIATIONS, BONFERRONI AND LORENZ CURVES 


Deviations from the mean and median help in giving a sense of the amount of spread 
in a population. The mean deviation about the mean and mean deviation about the 
median of the KLP distribution are given by 


oO 


6 B?w; 
D(u) = 2pF epr(u) — 24+ >) (B+ ram Tq L(8,0(k + 1), a,b, 1, u) 








j,k=0 
and 
i = 03? Ww; % 
D(M) aac! ae 2 (B + 1)[e%4+) Sy ee ree), 
J; i 


Consequently, Lorenz and Bonferroni curves are given by 


Lg) = KDE, and BC (a)) = Masel) 


or 
1 q 1 q 
L(p) T = iene (t)dt, and B(p) = =f Lies (t)dt, 


respectively, where q = F71 (p). 


KLP 
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6. MAXIMUM LIKELIHOOD ESTIMATION 


Let 71,--- ,£n be a random sample from the KLP distribution. The log-likelihood 
function is given by 


L = nlog(a) + nlog(b) + nlog(@) + 2nlog(G) — nlog(6 + 1) 
+ X log(l+z:) -8X r+ A+ (a-1) X log(e* -— 1) 
i=1 i=0 i=1 =j 


+ (b-1) > log (e — 1)° — (e™ — nel — nablog(e? — 1). 





The elements of the score vector are given by 





on = + oe — 1) — nblog(e?® — 1) 
(e? — 1)*log(e® — 1) — (e™ — 1) log(e* — 1) 
+ O-DY (e — 1) — (e — 1) 
OL n 


E AS On = a Mi a} Oo 
Et + a C 1)* —(e 1) | na log(e’ — 1), 





ði Ai 
OL “. 36 


3B = => Saag -Jard Ge (a — 1) DD en 
^M 1) 10d; 


= a(b — 1) > (e? 5 1)4 = (er 2 











and 
OL n 2 tV xije™ nabe? 
æ ə > V (xi) ; 








+ 








Orn Be p WOON 1 
g7” G5) aa 


and 





The maximum likelihood estimates, A of A = (8,0,a,b) are obtained by solving the 


: : l palgal _ at . . 
nonlinear equations 35 = 9356 = 954 = 9, and 2 = 0. These equations are not in 
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closed form and the values of the parameters 6,0,a and b must be found by using 
iterative methods. 

We maximize the likelihood function using NLmixed procedure in SAS as well as 
the function nlm in R ([32]). These functions were applied and executed for a wide 
range of initial values. This process often results or leads to more than one maximum, 
however, in these cases, we take the MLEs corresponding to the largest value of the 
maxima. In a few cases, no maximum was identified for the selected initial values. In 
these cases, a new initial value was tried in order to obtain a maximum. 

The issues of existence and uniqueness of the MLEs are of theoretical interest and 
have been studied by several authors for different distributions including [28], [27], [34], 
and [33]. At this point we are not able to address the theoretical aspects (existence, 
uniqueness) of the MLE of the parameters of the KLP distribution. 

Note that the KLP density fxzp(-; A) has second derivatives with respect to the 
parameters, so that Fisher information matrix (FIM), I,;(A) can be expressed as 

O log(fre(X;A4))] . . 
1(A) = a | 35.05, jagaaaaa 
Elements of the FIM can be numerically obtained by MATLAB or MAPLE software. 
The total FIM I,,(A) can be approximated by 
2 o? log L 
(6.1) J,(A) = | 95:05; 


For real data, the matrix given in equation (6.1) is obtained after the convergence 
of the Newton-Raphson procedure in MATLAB or R software. Let A = (8,0,â,b) 


be the maximum likelihood estimate of A = (6,6,a,b). Under the usual regularity 
conditions and that the parameters are in the interior of the parameter space, but not 


on the boundary, we have: /n(A—A) Er N,(0,1-'(A)), where I(A) is the expected 
Fisher information matrix. The asymptotic behavior is still valid if I(A) is replaced 





4,9 =1,2,3,4. 





al 4x4 


by the observed information matrix evaluated at A, that is J(A). The multivariate 
normal distribution with mean vector 0 = (0,0,0,0)7 and covariance matrix I~1(A) 
can be used to construct confidence intervals for the model parameters. That is, the 
approximate 100(1 — 7)% two-sided confidence intervals for 3,0,a and b are given by 


B E 27/2 V 1,3(A), ô+ Zn/2 V Tip (A), a+ Zy/2 V RHA); 


and 6 + Zn/2 Iz (Â), respectively, where T53(A), Tog (A), 13 (A) and I} (Â) are 
diagonal elements of I71(A) = (nI(A))~! and Z,/2 is the upper (7/2) percentile of a 
standard normal distribution. 

We can use the likelihood ratio (LR) test to compare the fit of the KLP distribution 
with its sub-models for a given data set. For example, to test a = b = 1, the LR statistic 
is w* = 2[In(L(G, ô, G, b)) — In(L (3,0, 1,1))], where 8, Ô, @ and 6 are the unrestricted 
estimates, and § and @ are the restricted estimates. The LR test rejects the null 
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hypothesis if w* > x?, where x? denote the upper 100e% point of the x? distribution 
with 2 degrees of freedom. 


6.1. Monte Carlo Simulation Study. In this sub-section, we study the perfor- 
mance of the maximum likelihood method for estimating the KLP model parame- 
ters by conducting simulations for different sample sizes and different parameter val- 
ues. Equation in (2.9) was used to generate random data from the KLP distribu- 
tion. The simulation study was repeated N = 1,000 times each with samples of size 
n = 100, 200, 400, 800, 1000 and parameter values J : 8 = 0.5,0 = 0.4,a = 0.3,b = 0.5 
and JI : B = 2.0,0 = 2.0,a = 0.5,b = 0.5. Four quantities were computed in this 
simulation study: 


(a) Average bias of the MLE Ô of the parameter 0 = 6,6,a,D: 
La 
a SÓ- 9). 
i=1 


(b) Root mean squared error (RMSE) of the MLE Ô of the parameter V = 8,0,a,b: 





i=1 
(c) Coverage probability (CP) of 95% confidence intervals of the parameter 0 = 
B,0,a,b, i.e., the percentage of intervals that contain the true value of the 
parameter V. 
(d) Average width (AW) of 95% confidence intervals of the parameter 0 = 6,6, a,b. 
Table 6.1 presents the Average Bias, RMSE, CP and AW values of the parameters 
G8,0,a and b for different sample sizes. According to the results, it can be concluded 
that as the sample size n increases, the RMSEs decay toward zero. We also observe 
that for all the parameters, the biases decrease as the sample size n increases. Also, 
the results show that the coverage probabilities of the confidence intervals are quite 
close to the nominal level of 95% and that the average confidence widths decrease as 
the sample size increases. 
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TABLE 6.1. Monte Carlo Simulation Results: Average Bias, RMSE, CP 
and AW 




















I II 
Parameter n Average Bias RMSE CP AW Average Bias RMSE CP AW 
B 100 -0.00501 0.24462 0.9510 1.51624 -0.10803 0.87209 0.9410 5.20699 
200 -0.00295 0.23610 0.9490 1.13275 -0.03140 0.76234 0.9530 4.25096 
400 -0.00043 0.20765 0.9480 0.91630 -0.02458 0.64276 0.9580 3.24892 
800 -0.00027 0.18282 0.9450 0.70245 -0.01504 0.49704 0.9550 2.38966 
1000 -0.00019 0.16592 0.9550 0.66421 -0.00511 0.45924 0.9580 2.08619 
0 100 0.88197 2.34717 0.9610 9.92480 0.84490 1.69587 0.9520 7.14813 
200 0.65401 1.90481 0.9520 7.38017 0.40583 1.08716 0.9480 4.51975 
400 0.61498 1.35036 0.9490 5.00112 0.18632 0.74348 0.9450 3.01703 
800 0.47330 0.96665 0.9460 3.40629 0.07429 0.49949 0.9520 2.02828 
1000 0.38674 0.79037 0.9510 2.87423 0.05577 0.45991 0.9490 1.80380 
a 100 0.00487 0.05577 0.9480 0.29526 -0.00084 0.13179 0.9370 0.59255 
200 0.00189 0.04795 0.9490 0.21880 -0.00056 0.10147 0.9420 0.43481 
400 0.00046 0.03968 0.9530 0.16529 -0.00019 0.08019 0.9450 0.31945 
800 0.00037 0.03462 0.9610 0.12181 -0.00013 0.05631 0.9480 0.22884 
1000 0.00016 0.02995 0.9530 0.11095 -0.00103 0.05113 0.9530 0.20353 
b 100 0.89313 1.65094 0.9450 10.26548 0.77219 1.75168 0.9640 10.59793 
200 0.88861 1.54050 0.9490 7.79171 0.48749 1.11044 0.9540 5.20265 
400 0.60264 1.13678 0.9510 4.38831 0.29611 0.72701 0.9480 2.79117 
800 0.42173 0.75993 0.9420 2.50884 0.13845 0.34273 0.9560 1.21017 
1000 0.31041 0.56817 0.9480 1.84970 0.12892 0.30152 0.9510 1.03431 





7. APPLICATIONS 


In this section, the KLP distribution is applied to real data sets in order to illus- 
trate the usefulness and applicability of the model. We fit the density functions of the 
Kumaraswamy Lindley-Poisson (KLP), Lindley-Poisson (LP) and Lindley (L) distri- 
butions. For comparison purposes, we also fit the beta exponentiated Lindley (BEL) 
distribution [22] which is also a 4 parameter model comparable to the KLP distribution. 
The pdf of the BEL distribution is given by 

forh; 3,0, a,b) BE O (1+ aje 
T, P, U, Q, = + x )e 
cor B(a,b)(8 +1) 
Bat 








x 1— {1+ 


6g b-1 
px e Bt 


6+1 


for x > 0,8 > 0,0 > 0,a > 0,b > 0. Estimates of the parameters of the distributions, 


A 


standard errors (in parentheses), Akaike Information Criterion (AJC = 2p — 2 log(L)), 
Consistent Akaike Information Criterion (AICC = AIC + 2+0), Bayesian Infor- 


n—p-1 
mation Criterion (BIC = plog(n) — 2log(L)), where Ê = L(A) is the value of the 
likelihood function evaluated at the parameter estimates, n is the number of observa- 
tions, and p is the number of estimated parameters are obtained. 

The first data set consists of 119 observations on fracture toughness of Alumina 
(Al,O3)(in the units of MPa mt). This data was studied by Nadarajah and Kotz 
[19]. The data is available at and taken from the following web-site: 
http://www.ceramics.nist.gov/.srd/summary/ftmain.htm. 


x 1-—|J1l1-—{1+ 
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The second data set gives failure times of a sample of n = 101 aluminum specimens 
of type 6061-T6 obtained by Birnbaum and Saunders [4]. These specimens were cut 
parallel to the direction of rolling and oscillating at 18 cycles per seconds and they were 
exposed to a pressure with maximum stress of 31,000 pounds per square inch (psi). 
The specimens were tested until failure. 

The third data set that is fitted to the KLP distribution consists of breaking stress 
of carbon fibers which was analyzed by Bader and Priest [2]. The data represent the 
tensile strength, measured in GPa, of 69 carbon fibers tested under tension at gauge 
lengths of 20 mm. 

The fourth data set set consists of 63 observations of the strengths of 1.5 cm glass 
fibres, originally obtained by workers at the UK National Physical Laboratory. The 
data was also studied by Smith and Naylor [31]. 

Table 7.1 displays results obtained from analyzing the Silicone Nitride data studied 
by Nadarajah and Kotz [19]. Estimates of the parameters of KLP distribution (s- 
tandard error in parentheses), Akaike Information Criterion (AIC), Consistent Akaike 
Information Criterion (AICC), Bayesian Information Criterion (BIC), Kolmogoroy- 
Smirnov (K-S) statistic and its p-value are given in the Table 7.1. Plots of the fitted 
densities and histogram, and observed probability versus predicted probability for the 
Silicone Nitride data are given in Figures 7.1 and 7.2, respectively. For the probability 


Sa 20475 
plot, we plotted Fxip(x(;); 6,9, a,b) against eer j=1, 2, --+,n, where £j are 
n 


the ordered values of the observed data. The measures of closeness are given by the 
sum of squares 
n j —0.375\]* 
SS = J F >) — | ——__] ; 
a | xoP(%)) ( n +0.25 


TABLE 7.1. Estimates of Models for Silicone Nitride Data 





Estimates Statistics 
Model B 0 a b -2logZL AIC AICC BIC S$ K-S p-value 
KLP(6,6,a,6) 0.3310 4.3514 2.2709 108.36 336.9 344.9 345.3 356.0 0.0651 0.063963 0.7150 
(0.2214) (5.9017) (2.2192) (339.92) 





LP(8,6,1,1) 11076 22.9292 1 1 368.3 3723 3724 377.8 0.5760 0.142930 0.01547 
(0.0623) (4.7641) 

L(3,-,1,1) 0.3967- 1 1 534.8 536.8 536.8 539.6 4.9662 0.335270 4.812x10 
0263) 





BEL(8,,a,b) 0.2120 6.0217 0.6111 431.95 337.0 345.0 345.4 356.2 0.0784 0.070035 0.6037 
(0.1067) (4.6244) (0.4978) (49.4370) 














The LR test statistics of the hypotheses Hp : LP against Ha : KLP, Ho: L against 
H, : KLP, and Ho: L against H, :LP are 27.4 (p-value=1.12x10~® < 0.001), 197.9 (p- 
value=1.20x107/? < 0.001) and 166.5 (p-value=4.30 x10738 < 0.001). We can therefore 
conclude that there is a significant difference between KLP and LP distributions, KLP 
and L distributions as well as between LP and L distributions. The values of the 
statistics AIC, AICC and BIC are very close for the KLP and the BEL distributions. 
However, the KLP distribution has the smallest K-S statistic and the largest p-value 


KUMARASWAMY LINDLEY-POISSON DISTRIBUTION 19 


indicating that KLP distribution provides a better fit among the four models when 
fitting the Silicone Nitride data. The KLP yields the smallest value for Sum of Squares 
(SS) among all the models. 


TABLE 7.2. Estimates of Models for Aluminum Specimens Data 





Estimates Statistics 
Model B 6 a b -2logL AIC AICC BIC SS K-S p-value 
KLP(3,6,a,b) 0.0313 6.1948 3.9159 4.3329 904.4 912.4 912.8 922.8 0.0413 0.061896 0.8382 
(0.0062) (23.8570) (14.5997) (2.4939) 





LP(8,0,1,1) 0.0532 96.3691 1 1 914.7 918.7 918.8 923.9 0.1759 0.100170 0.2682 
(0.0033) (32.6069) 

L(B,-,1,1) 0.0148 - 1 1 1106.3 1108.3 1108.3 1110.9 4.8973 0.400330 2.398x10-! 
(0.0011) 


BEL(,0,a,b) 0.0318 31.6377 0.6532 42752 9049 912.9 9123 923.3 0.0465 0.067242 0.7565 
(0.0068) (98.0430) (2.2797) (6.0065) 














Table 7.2 gives the estimates of the model parameters and the statistics AIC, AICC, 
BIC, K-S statistic and its p-value for the Aluminum specimen data. The LR test 
statistics of the hypotheses Hp : LP against H, : KLP, Hp : L against Ha : KLP, and Ho : 
L against Ha :LP are 10.3 (p-value=5.8x 107? < 0.01), 201.9 (p-value=1.64x10~* < 
0.001) and 191.6 (p-value=1.42 x10743 < 0.001). We can therefore conclude that there 
is a significant difference between KLP and LP distributions, KLP and L distributions 
as well as between LP and L distributions. The values of the statistics AIC, AICC 
and BIC are very close for the KLP and the BEL distributions. However, the KLP 
distribution has the smallest K-S statistic and the largest p-value indicating that KLP 
distribution provides a better fit among the four models when fitting the Aluminum 
Specimens Data. The KLP yields the smallest value for Sum of Squares (SS) among 
all the models. Plots of the fitted densities and histogram, observed probability versus 
predicted probability for the Aluminum specimen data are given in Figures 7.3 and 
7.4. The plot of the fitted pdf support the conclusion based on Table 6.2. The figures 
suggest that both KLP and BEL distributions captures the middle part of the data, 
as well as the tails better than the fitted sub-models. 

Table 7.3 shows the results obtained from analyzing the carbon fibers data of Bader 
and Priest [2]. 


TABLE 7.3. Estimates of Models for Carbon Fibers Data 





Estimates Statistics 
Model b 0 a b -2log L AIC AICC BIC SS K-S p-value 
KLP(6,0,a,b) 0.7230 0.0348 8.3910 24.7633 97.7 105.7 106.3 114.7 0.0136 0.038759 0.9999 
(6.2713) (37.1709) (33.3868) (429.21) 





LP(8,0,1,1) 23266 67.1553 1 1 107.8 111.8 1120 116.3 0.1332 0.088736 0.6489 
(0.1808) (24.3584) 

L(8,-,1,1) 0.6545 : 1 1 -238.4 2404 2404 242.6 3.4678 0.401130 4.547 10-29 
(0.05803) 


BEL(8,0,a,b) 1.3696 48,7520 0.3328 7.2260 97.9 105.9 106.5 114.8 0.0161 0.041807 0.9997 
(0.4863) (99.3139) (0.5132) (2.0853) 
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The LR test statistic of the hypothesis Ho : LP against Ha : KLP, Hp : La 
gainst H, : KLP, and Hp : L against H, :LP are 10.1 (p-value=6.4x 107%, 140.7 (p- 
value=2.67 x 10~*° < 0.001) and 130.6 (p-value=3.03 x 10~*° < 0.001). We can therefore 
conclude that there is a significant difference between KLP and LP distributions, KLP 
and L distributions as well as between LP and L distributions. The values of the 
statistics AIC, AICC and BIC are very close for the KLP and the BEL distributions. 
However, the KLP distribution has the smallest K-S statistic and the largest p-value 
indicatiing that KLP distribution provides a better fit among the four models when 
fitting the Carbon Fiber data. The KLP yields the smallest value for Sum of Squares 
(SS) among all the models. 

Plots of the fitted densities and histogram, observed probability versus predicted 
probability for the Carbon Fiber data are given in Figures 7.5 and 7.6, respectively. 
The plot of the fitted pdf support the conclusion based on Table 7.4. The figures 
suggest that both KLP and BEL distributions capture the middle part of the data, as 
well as the tails better than the fitted sub-models. 


TABLE 7.4. Estimates of Models for Glass Fibers Data 





Estimates Statistics 
Model p 6 a b —2logL AIC AICC BIC S$ K-S p-value 
KLP(6,0,a,b) 0.7334 11.8621 1.2658 2083.13 28.2 36.2 36.9 44.7 0.1529 0.13306 0.2146 
(0.3463) (22.4170) (2.3915) (9883.45) 





LP(8,0,1,1) 3.0456 29.3432 1 1 59.3 63.3 63.5 67.6 0.7157 0.21887 0.004783 
(0.2338) (8.3871) 

L(8,—,1,1) 0.9961- 1 1 162.6 164.6 164.6 166.7 3.3017 0.38642 1.349x10-8 
(0.0948) 


BEL(8,,a,b) 0.4705 8.0340 0.5670 8989.04 29.2 37.2 37.8 45.7 0.1951 0.14755 0.1287 
(0.1995) (4.2325) (0.3193) (124.49) 





Plots of the fitted densities and histogram and probability plots for the glass fibres 
data from Smith and Naylor [31] are given in Figures 7.7 and 7.8 respectively. The 
LR test statistic for the test of the hypotheses Ho : LP against Ha : KLP, Ho : L 
against H, : KLP, and Hp: L against Ha :LP are 31.1 (p-value=1.765x10~", 134.4 (p- 
value=6.092x 107° < 0.0001) and 103.3 (p-value=2.88 x 10~*4 < 0.0001), respectively. 
We can therefore conclude that there is a significant difference between KLP and LP 
distributions, KLP and L distributions as well as between LP and L distributions. 

The values of AIC, AICC and BIC shows that the KLP distribution is a better 
model and the SS value is comparatively smaller than the corresponding values for 
the LP and L distributions. The values of these statistics (AIC, AICC, BIC) for 
the KLP distribution are very competitive when compared to those of the non-nested 
BEL distribution. However, the KLP distribution has the smallest K-S statistic and 
the largest p-value indicating that KLP distribution provides a better fit among the 
four models when fitting the Glass Fiber data. The plot of the fitted pdf support 
the conclusion based on Table 7.4. The figures suggest that both KLP and BEL 
distributions captures the middle part of the data, as well as the tails better than the 
fitted sub-models. 
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FIGURE 7.1. Histogram and Fitted Density for Silicone Nitride Data 
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FIGURE 7.2. Probability Plots for Silicone Nitride Data 
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FIGURE 7.3. Histogram and Fitted Density for Aluminum Specimens Data 
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FIGURE 7.4. Probability Plots for Aluminum Specimens Data 
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FIGURE 7.5. Histogram and Fitted Density for Carbon Fibers Data 
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FIGURE 7.6. Probability Plots for Carbon Fibers Data 
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FIGURE 7.7. Histogram and Fitted Density for Glass Fibers Data 
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FIGURE 7.8. Probability Plots for Glass Fibers Data 


8. CONCLUDING REMARKS 


We have proposed and presented a new class of lifetime distributions called the Ku- 
maraswamy Lindley-Poisson distribution. This class of distributions has applications 
in lifetime data analysis. The KLP distribution has the LP, exponentiated Lindley 
Poisson and Lindley distributions as special cases. The density of the new distribution 
can be expressed as a linear combination of Lindley Poisson distributions. The KLP 
distribution has a hazard function that displays flexible behavior. Moments, mean de- 
viations, distribution of order statistics and Rényi entropy were derived. The method 
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of maximum likelihood was used to estimate the model parameters. Finally, KLP dis- 
tribution is fitted to real data sets in order to illustrate the applicability and usefulness 
of the distribution. 
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